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Abstract. - We present a field-theoretic renormalization group (RG) analysis of a single 
flexible, screened polyelectrolyte chain (a Debye-Hiickel chain) in a polar solvent. We point 
out that the Debye-Hiickel chain may be mapped onto a local field theory which has the same 
fixed point as a generalised n ^ 1 Potts model. Systematic analysis of the field theory shows 
that the system is one with two interplaying length-scales requiring the calculation of scaling 
functions as well as exponents to fully describe its physical behaviour. To illustrate this, we 
solve the RG equation and explicitly calculate the exponents and the mean end-to-end length 
of the chain. 



Polyelectrolytes are polymers with ionizable monomers which in polar solvents, such as water, 
dissociate into charged polymers and small 'counter-ions' of opposite sign. They are of 
widespread importance with many applications in physics, biology and chemistry y. Their 
solublility in polar solvents means that they also have many industrial applications. Typical 
examples are DNA and sulphonated polystyrene. 

The behaviour of polyelectrolytes is quite poorly understood, chiefly because it has been 
difficult to deal with their long-range coulomb interactions theoretically. In addition, counter- 
ions, complex-formation and salt and their differing length-scales make comparison between 
theory and experiment difficult because of questions of validity of the theories in different 
regimes. Neutral polymers on the other hand, have been extremely succesfully described by 
scaling ideas and renormalization group theories ||^. The separation of the length-scales in 
their physics leads to universality. The main properties of the system are independent of 
microscopic details which only change prefactors in the physical quantities and make analysis 
of experimental results, at least in their limiting behaviour, reasonably straightforward. It 
would be particularly useful if the behaviour of polyelectrolytes could also be put in a similar 
context. The resulting theory would obviously be more complicated reflecting the increased 
complexity of the system. Nonetheless a 'universal' theory of polylectrolytes would be a real 
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step forward in the understanding of polymer solutions. This is what we attempt to do in 
this paper. In this spirit we consider the simplest possible system, a single eharged chain in 
solution. 

In solution the counter-ions in the vicinity of the macro-ion(polymer) 'screen' or reduce 
the range of the coulombic interaction. The simplest model of a screened polyeleetrolyte 
which we use as the starting point for our analysis is the Wiener chain with a Debye-Hiickel 
(DH) or Yukawa potential interaction between charged monomers. The chain is described by 
a vector r(s) parametrised by the arc-length s and the partition function is obtained from 
Z = J 2?r(s) exp {—(3Ti.DH [r(s)]} with the Hamiltonian of the system given in d dimensions by 

l^'^DH t ds {^^ +b r r dsds'V{v{s) - r(s')), (1) 

Jo Jo 



2e Jo \ds 



where i is the ^ Kuhn- segment length of the chain, L — £N the length of the chain, where 
N is the number of segments on the chain. The Kuhn-length here refers to the distance 
between charged monomers (see inset Fig.l). 1/(3 = ksT and b cx where 
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is the Bjerrum length which measures the strength of the bare coulomb interaction. In water, 
Xb = 7.14A. As usual eeo is the dieleetic constant of the solvent, ks is the Boltzmann constant 
and T the temperature. The potential between monomers, V{r) is a solution of the linearised 
Poisson-Boltzmann equation (—A -I- K^)V^(r) — which in three dimensions is given by the 
usual V{r) = e^''l''l/|r|. The screening length or range of the interaction is given by 1/n and 
K oc •'/Xb is a function of the density of screening ions and the dielectric properties of the 
solution. The DH model js] is thought to be valid only when the polyeleetrolyte is weakly 
charged (see inset Fig.l) and in the presence of salt. There are no large fluctuations in the 
counter-ion/salt density so for the Coulomb interaction the DH model may be considered to 
be ^ mean- field. 

Recent extensive simulations of flexible DH chains Q give dramatically different results from 
existing theories |^. We compare our results with these simulations because experiments 
are always performed at non-zero concentrations where many-chain effects begin to play a part 
and cannot be used to check the validity of a single chain model The quantity that we 
want to calculate is the average end-to-end length (i?^) = {{r{L) — r(0))^). 

Returning to equation (Q), we apply a generalisation of the well known de Gennes trick 
to map the polyeleetrolyte to an n — > field theory [§ with a non-local (long-ranged) 
interaction H, ^ . The propagator is defined in the usual way as the Laplace transform (LT) of 
the 2-point function G'(R) = /[X>0]0(R)<?!)(O)e-^('^) = dCe-*^ {S{rlc) ~ r(0) - R)) where 
C = 2Le/d = 2Ne'^/d with the action given by 5((/.) = [^{^(t)? + X h' (t> {i^)V{r- 

r')(fP{r') where 2g — {2£/d)'^b and 0(r) are n ^ component fields with 0{n) symmetry. We 
have used the notation = J d'^r. We perform a Hubbard-Stratonovich transformation to 
remove the DH potential and generate a new scalar field V'(r) so that 

G(R) = j -^e^^-G{q)=M I mivmmWe-'^'^^^^ (2) 



with S{4>, V) 



(3) 



where 4"'' is the polyeleetrolyte fugacity, ij"-' = and A/" is a normalisation factor. The 
coupling is given by wo/3! = i^i/g) if ip is real(imaginary). ip{r) is a fluctuating Coulomb 
field which is coupled to the polymer density. It has fluctuations on a length-scale given by 
1/k. Since there are two independent fields, we have two sets of exponents i^^, i^,/,; ^y^, ry^. To 
distinguish between the two we use a superscript ip for the Coulomb terms and a superscript 
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(f) for the polymer terms. 

There is a rich literature on the use of the RG and field theory to calculate the properties 
of neutral polymers. On the other hand, these techniques have not been applied successfully 
to charged polymers, polyelectrolytes. This was due to a misconception that they could not 
be usefully applied to such systems in the physical dimension, d — 3. In the literature on field 
theory and renormalization group theory of polyelectrolytes |p, , work has always been 
done on unscreened chains. This has led to an unfortunate obscuring of the physics which we 
believe we have finally clarified. The exponents can be calculated to all orders in e = 6 — d 
(the upper critical dimension d^ ~ 6 can be inferred from the Lagrangian (||)) and is given by 

— 2/ [d — 2) H]. It has always been assumed that the end-to-end length is given by 

and since = 1 in d = 4, and the chain could not be stretched out more than a rod , 
{E?) PN"^, it was said that due to overstretching of the Gaussian chain the theory breaks 
down for d < A. It was therefore thought that field theoretic methods could not be used on 
charged chains in the physical dimension d = 3. We will show that this analysis is incomplete. 

We have a number of new results in our letter. Firstly, we point out a connection between 
a Potts model used to study percolation and the sol-gel transition and a single polyelectrolyte 
chain |TT| . Secondly, we include the two important length-scales which are necessary for 
the analysis of the system. They are the screening-length 1/k and the gaussian chain size 
S = (Rq)^^^ ~ iN'^/^. A wide range of systems can be studied beginning where the screening 
is very high (SAW), K^t'^N >> 1 (the polymer field </> is more critical than coulomb field 
ip) and ending where it is very low (unscreened, rod- like), n'^i'^N << 1 (ip is more critical 
than (/)). This is therefore a problem of cross-over (l2j. The previous work has focussed, in 
our opinion mistakenly, on the unscreened case by setting n — 0. The exponent is actually 
calculated for the intermediate regime ~ 1/{£^N) — > and does not apply to the unscreened 
(rod) regime where k and Hi^PN — > 0. We also find that quantities like the end-to-end 
distance depend on negative powers of k. This implies that the operator k^^"^ in eqn. (|^) 
is dangerously irrelevant . By setting k = at the start of the calculation, one loses sight 
of such subtleties. 

We calculate in agreement with earlier work ^, the exponents 

e 1 2 

where e — Q — d. The exponents can also be obtained from the percolation theory but at a 
different fixed point from the usual sol-gel transition This leads naturally to our final 

conclusion: these exponents correspond to a bi-critical fixed point (both (/), -ip are critical) 
where the screening length and the chain size are comparable and at this fixed point the 
end-to-end length we calculate as 

where R/t — k and f{x) is a function singular as x —> 0. We explicitly calculate f{x) = 
^(-2^^+i)^(-^) where g{x) = (1 - 3elog(a;)/(4a;))[l - e(37/(4a;) - 5/24)][l - e(3/(4x) - 1/2 + 
a;/8) exp (x)i?i (a;)]-|-0(e2). The exponential integral i?i (2:) is defined by i?i(z) = dtexp{~zt)/t 
and 7 = 0.57721 is Euler's constant. The end-to-end length is given by 

(Recall that at this fixed point ~ l/{i'^N) 0). Earlier work ^ had always implicitly 
assumed that f{x) was a well behaved function in the limit x [k 0). We find 
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2i'[ = 2 — l/{2i>^) determines the chain size and not 2v^. Also, because of g the exponent 
will vary slightly with n [Q. There is consequently no over-stretching (at this fixed-point, 
u*^ — — 27e/4) for d < 4 and equation (Q) should be replaced by equation (Q). It is important 
to note that v[ — 7/8 = 0.875 (intermediate between the SAW, v — 0.588 and rod, v — 1.0) in 
three dimensions. It is also clear from our analysis that the lower critical dimension is d = 2 
where the theory breaks down. 

In order to compare with earlier approaches we may 'fit' our results from equation (|^) to a 
worm-like chain model jl3|. The electrostatic persistence length, le, may be estimated from 
€e(l — e~^^/^=) cx R'^^g[R'^ N'^) with ^ — —3/2 in d = 3 showing no unique power law behaviour. 
This is in agreement with the simulations of Micka and Kremer The rod- like behaviour 
in the long-range, low salt regime 1/k^ >> f'N and the self-avoiding walk behaviour for the 
short-range, high salt regime, i'^N >> [|i^ can be obtained using different methods 

Now we perform a RG analysis of our model with Lagrangian given by equation (^). Our 
goal is to calculate 

= -T-\d,2G{q)\,=o)/'rc\G{q)\,^o) - D{u)F{£,h) (8) 

where D{u) is a non universal constant, is the inverse LT with respect to C and G{q) = 
(7(2,0) is the FT of the renormalised chain propagator. The bare parameters are expressed 
in terms of their renormalized values by uq = ^'^^'^ SJ^ ZuU, t^^""* = IJ^Z^^ [i2-^-!/'.0^ + ^i-^-i/..^/.^] 

and 4°'' = fJ'^Z^^ [^2.^0, 02 + tjZ^ ,^2] where is an external length-scale and Sd = fld/i^TrY 
with fid the angular part of a d-dimensional integral. We denote an N point polymer M point 
coulomb vertex function as r*^^'*^^ while a composite N point polymer M point coulomb 
vertex function with L, ifP' and ip'^ insertions is generally given by r(^>^^*^^^). The 
wave-function rcnormalization factors are calculated by expressing the renormalized vertex 
functions in terms of the bare ones as r^'*^'' = z^^^ Z^^^^T^^'^^\ The Z-factors are evaluated 
using the vertex functions which are calculated using a dimensional regularisation scheme with 
minimal subtraction of poles jl^. We begin with the Coulomb vertex F^^'^^ which due to the 
n —f limit gives the result = 1 because all closed subdiagrams composed of these polymer 
lines give a factor n. From the polymer vertex F^^'*^^ we get Z^. For the determination of 
Zu we need the leading singular 1/e contribution of F^^'^^ from which follows the /3-function 
the zero of which yields the infra-red stable fixed point w*^ — — 27e/4. An evaluation of 
r(o,0;2,i) yigi(js Z^^^2 = 1 and r^"'^'^'"' gives Z^^^2 = using similar arguments as those for 
Z^. The composite polymer vertex function r^^'^'^'^^^ gives Z^^^2 and F^^'"-"'^^ to 1-loop yields 

Zcji^^^ = ^0,02 — 1. 

To see the scaling form of (i?^), we need to solve the RG equation for Q{q) = T^^[G^^'^\q)] 
which is given by 

{fid^ + (iudu - 2tidt, - ^2Cdc + 70 - ^2 - ^ihC) g{q/fi; u, £, h) = 0. (9) 

where f3u = {du/d\nfi), 7^ = /3„9„lnZ^, 7^ = /3„9„lnZ0, 7^2 = In ^^^^2, ??i = 

— ^0/^0 02/3„9„[Z0 ,^2 /Z^Jand — —2— 7^2-1-70 and we have used the notation dxA = dA/dx. 
The solution is readily obtained using the method of characteristics [|4| giving 

F{CM) = /:'"*Fi(l,ti/:2''*,?/*) = l/iiF2(£ty^'"*\l,u*) (10) 

showing non-trivial scaling (t]^ — — 2 = 7^2(7/*) — 70(7/*)). This means that the 

scaling function must be explicitly calculated to extract information about the size of the 
chain. Great care must be taken in the interpretation of the inverse Laplace transform. As we 
are dealing with critical fields, the masses ti and t2 must go to zero. We get non-trivial results 
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Fig. 1. - We plot Micka and Kremer's 3d simulation data of flexible DH chains for different values 
of I, N and k. ^{i?2)^2-i/,.^y^ kV"* jv with = 2/(d - 2) evaluated at d = 3 and we get all 
the points to collapse on one curve. Inset: weakly charged polyelectrolyte and Kuhn length t. 



only when both length-scales are large and of the same order otherwise the larger length-scale 
dominates the physics and the smaller becomes irrelevant. This means that the ratio of the 
two length-scales c = must be finite. This gives the non-trivial fixed-point u* ^ 0. ft 
is reassuring to note that the explicit calculation, (|^), agrees exactly with the general form of 
the scaling function in equation (p^). 

As is usual in RG analysis one must check for stability of the fixed point as the quartic terms 
become relevant as we cross the dimension d — A. This says something about the validity of 
the DH model (which assumes these operators are irrelevant) in the physical dimension d = 3. 
Any operators which make the fixed-point unstable are corrections to DH theory. The -0^ 
and ^'^(j)^ and (0^)^ are allowed by symmetry and become relevant for d < 4 |9j. The excluded 
volume term {(j)^)'^ is only relevant in the high salt regime. The quartic corrections mean that 
the DH model is invalid below d = 4 in the limit k ^ because the critical fluctuations of 
the Coulomb field for d < 4 mean that it cannot be described by a Gaussian model (no terms 
higher than quadratic in t/;) [ p^ . Since we know that DH theory is in a sense a 'mean-field' 
approximation, we should expect it to be wrong in low dimensions. Therefore is valid for the 
DH model for d < 4 but not for the complete theory (DH + corrections). The approximation 
that the fluctuations in the density of the screening ions are small breaks down. This does 
not really have anything to do with > 1, d < 4 which is a result of a field theory of a DH 
chain. 

Flexible polyelectrolytes have been simulated by Stevens and Kremer ||l^ using explicit 
counter-ions and Micka and Kremer Q using a DH potential. Stevens and Kremer found 
exponents vi ~ 0.94 < 1 in good agreement with our picture. Micka and Kremer Q observed 
no unique power law as k and L were varied. In Fig.l we plot their end-to-end data using our 
explicit scaling form in equation(^ and get very good agreement. More detailed comparison 
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of this and other quantities with their numerical data wiU be presented elsewhere [|17[ . 

In conclusion, we have systematically calculated the properties of a single screened poly- 
electrolyte in the Debye-Hiickel approximation (a well-defined problem) and found non-trivial 
scaling behaviour in good agreement with simulation. The DH model is therefore controlled by 
three fixed points governed by c = PNk"^. They are c — > oo (SAW), c ~ 1 the new non-trivial 
behaviour which is the subject of this paper and c — > (rod). Our analysis confirms that the 
DH model breaks down if there are large fluctuations in the coulomb field, as would be observed 
in counter-ion condensation. We hope that this paper will stimulate the interest in a 20 year old 
problem in the literature, our conclusion being that the calculation of the exponents by Pfeuty 
et al and others |, |, was absolutely correct but their interpretation of those exponents 
and hence of the physics was incomplete. We believe that this paper therefore re-opens 
the area of the field theoretic RG description of charged polymers. In future publications 

we will deal with charged manifolds, directed polyelectrolytes, solutions and networks [ p^ . 
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